A simple approach for charge renormalization of highly charged macro-ions 
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We revisit the popular notion of effective or renormalized charge, which is a concept of central 
importance in the field of highly charged colloidal or polyelectrolyte solutions. Working at the level 
of a linear Debye-Hiickel like theory only, we propose a simple, efficient and versatile method to 
predict the saturated amount of charge renormalization, which is however a non-linear effect arising 
at strong electrostatic coupling. The results are successfully tested against the numerical solutions of 
Poisson-Boltzmann theory for polyions of various shapes (planar, cylindrical and spherical), both in 
the infinite dilution limit or in confined geometry, with or without added electrolyte. Our approach, 
accurate for monovalent micro-ions in solvents like water, is finally confronted against experimental 
results, namely the crystallization of charged colloidal suspensions and the osmotic coefficient of 
B-DNA solutions. 



Our present understanding of charged macro-ions sus- 
pensions is essentially based on the DLVO theory, named 
after Derjaguin, Landau, Verwey and Overbeek [1]. This 
approach relies on a Poisson-Boltzmann (PB) mean- field 
description of the micro- ions clouds. An important pre- 
diction of the theory is the effective interaction pair po- 
tential between two macro-ions (e.g. colloids) in the sol- 
vent which, within a linearization approximation, takes 
the well-known Yukawa or Debye-Hiickel (DH) form: 
v(r) ~ Z 2 exp(— nr)/r, where Z is the charge of the ob- 
ject and k denotes the inverse Debye screening length. 
However, this approach becomes inadequate to describe 
highly charged objects for which the electrostatic en- 
ergy of a micro-ion near the macro-ion surface largely 
exceeds ^T, the thermal energy, and the linearization 
of PB equations is a priori not justified. In this case how- 
ever, the electrostatic potential in exact [2] or mean-field 
[3-5] theories still takes a Debye-Hiickel form far from the 
charged bodies, provided that the source of the potential 
is renormalized [Z — > Z c g). The essential idea is that 
the micro-ions which suffer a high electrostatic coupling 
with the macro-ion accumulate in its immediate vicin- 
ity so that the decorated object (macro-ion plus captive 
counter-ions) may be considered as a single entity which 
carries an effective charge Z c g, much lower (in absolute 
value) than the structural one. Within PB theory, Z and 
Z c s coincide for low values of the structural charge, but 
Z e g eventually reaches a saturation value Z ssA indepen- 
dent of Z when the bare charge increases [5,6]. 

Of course, the difficulty remains to predict Z sat for a 
given suspension of macro- ions [4-7] . In the absence of 
any general analytical framework for the computation of 
the effective charge, this quantity is often considered as 
an adjustable parameter to fit experimental data [8,9]. 
In this Letter we show that a simple physical argument, 
at the level of the DH linearized description only, yields 



explicit (and in some favorable cases analytical) expres- 
sions for the effective charges at saturation, which com- 
pare well with both numerical solutions of non-linear PB 
theory and available experimental data. 

For simplicity, we start by considering a unique highly 
(positively) charged sphere immersed in a symmetric 1:1 
electrolyte of bulk ionic strength I = k 2 /(8tt£b), where 
£b = e 2 / (ATreksT) is the Bjerrum length (e being the 
dielectric constant of the solvent): £b — 7 A. for water at 
room temperature. Within PB theory, the dimensionless 
electrostatic potential $ = eV/ksT obeys the relation 



V $ = k sinh<F 



(1) 



Far from macro-ion (where it is understood that van- 
ishes), the solution $ pB of Eq. (1) also obeys the 
linearized Poisson-Boltzmann (LPB) equation V 2 $ = 
k 2 $, and therefore takes the Yukawa form $ LPB — 
Z c s£b exp[«:(a — r)]/[r(l + «a)], with a the radius of 
the sphere. Z e ff (in e units) is consequently defined 
here without ambiguity from the far field behaviour of 
<i>p B (see [4,10,11] for alternative definitions of effective 
charges). Accordingly, a "non linear" region may be 
defined (r £ [a,r*]), corresponding to $ pB larger than 
unity, where by definition of the cutoff r* , $ LPB {r* ) is 
of order 1. In the limit of large na, this "non-linear" 
region is however confined to the immediate vicinity of 
the macro-ion: r* ~ a. We consequently have the ef- 
fective boundary condition <!> LpB (a) ~ <& pB (r*) = C, 
where C is a constant of order 1, which yields immedi- 
ately Z e g = Ca(l + ko)/£b- This argument assumes that 
the bare charge Z is high enough to have <E> pB larger than 
unity close to the macro-ion, and therefore provides the 
saturation value of Z e g, denoted hereafter as Z sat [12]. 
We therefore easily obtain the non trivial dependence of 
this quantity upon physico-chemical parameters. 

This picture of a decorated macro-ion -where the 
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"bound" counterfoils renormalizing the charge appear 
to have an electrostatic energy eVo balancing the ther- 
mal energy fc^T- may be rationalized as follows. In the 
limit of large na, we perform an asymptotic matching 
of the non- linear PB planar solution (see [1]) to the lin- 
ear solution $ LPB in curved geometry. We obtain for 
high bare charges the same value of the contact po- 
tential ^ LPB (o) = 4 (of order 1 as expected) so that 
Z s &t = 4a(l + k<j)/£b [13]. Such a procedure provides 
by construction the correct large na (low curvature) be- 
haviour of Z sa t, but we will show below that it remains 
fairly accurate down to na of order 1. 

Generalizing this approach, we consequently obtain 
the leading curvature saturated effective charge from the 
following analysis. For an isolated macro-ion of arbi- 
trary shape in an electrolyte: a) find the electrostatic 
potential, 3> LPB , solution of the linearized PB equa- 
tions, supplemented by a fixed potential boundary con- 
dition: $ LpB (surface) = C, where C = 4 at leading or- 
der in curvature; b) deduce Z sat from Gauss theorem 
at the surface of the object. In the case of an infi- 
nite cylinder (radius a, bare lineic charge A), we obtain 
A sat = 2{Ka/ Ib)K\{ko) / Kq{ko), where Kq and K\ are 
the modified Bessel functions of orders and 1. 




/ca 



FIG. 1. Effective charge at saturation of an isolated spher- 
ical macro-ion (radius a) as a function of na. The continuous 
line is the analytical expression given in the text, while the 
dots are the results extracted from the far field behaviour of 
the non-linear PB potential. In the inset, the same results for 
the cylinder geometry are shown on a log-linear scale. 

In order to test the validity of our results, we have 
numerically solved the non- linear PB equation (1) for 
high Z values corresponding to the saturation regime, 
and computed the effective charge from the electrostatic 
potential at large distances (i.e. the value required to 
match $ LPB to the far field $ pB obtained numerically). 
Figure 1 compares the resulting PB effective charge to 
our expressions, for spherical and cylindrical macro-ions. 
The agreement becomes excellent at large na as it should, 
and in the case of cylinders, even holds down to very 



small na (0.01), a point which is not a priori expected. 
Finally, in the planar geometry our approach provides 
by construction the correct effective charge (compared 
to PB). 
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FIG. 2. Effective saturated charge of spherical macro-ions 
(radius a) as a function of volume fraction n for ko = 2.6. The 
continuous line shows the effective charge computed using the 
prescription, while the dots are the results of the non-linear 
PB theory, following Ref. [5]. Inset: no-salt situation. 

For spherical colloids, expressions reminiscent of that 
reported above can be found in the literature [4,14,15]. 
It seems however that the generality of the underlying 
method has been overlooked. In particular, our proce- 
dure may be extended to the finite concentration case, 
using the concept of Wigner Seitz (WS) cell [5] : the in- 
fluence of the other macro-ions is accounted for by con- 
fining the macro-ion into a cell, with global electroneu- 
trality. The size of the cell, Rws, is computed from the 
density of macro-ions, while its geometry is chosen as to 
mimic the spatial structure of the macro-ions in the so- 
lution. Suppose that the system is in equilibrium with a 
reservoir of salt defined in terms of its Debye length k _1 . 
We first linearize Eq. (1) around the value of potential 
at the cell boundary, <j>^ = &(R\ys), unknown at this 
point, which yields 

V 2 <5$ = nl ( 7o + <5$) (2) 

where (5$ = $ — $ Rl nl — K 2 cosh($ fl ) and 70 = 
y/l — (k/k*) 4 . Note that the relevant screening length 
k" 1 (always smaller than k~ 1 , a general feature for finite 
concentration) is not a parameter, and should be deter- 
mined at the end of the calculation. Equation (2) is sup- 
plemented by two boundary conditions: the consistency 
constraint [S^(Rws) = 0] and the global electroneutral- 
ity (which imposes a vanishing normal electric field at the 
WS boundary). To generalize the approach discussed in 
the limit of infinite dilution, we propose the following 
prescription (providing a third boundary condition) : the 
difference of potential between the macro- ion and the WS 
surface is <5$(a) = 4. Here again, the effective charge is 
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obtained from Gauss theorem at the macro-ion's surface. 

This generalized procedure is now explicited for a so- 
lution of spherical macro-ions with concentration p. The 
radius of the WS spherical cell is given as Rws = 
(47T/9/3) -1 / 3 . In this geometry, the (LPB) solution of 
Eq. (2) reads 



5$(r) = 70 



+ +/-- 



(3) 



where f± — {n±a ± I)/ (2k+) e~xjp(^FK*Rws)- Our pre- 
scription allows to compute k*, such that <5$(a) = 4. 
This equation is solved numerically for using a sim- 
ple numerical Newton procedure. The effective charge 
follows from the gradient of <5<f>(r) in Eq. (3) taken at 
r = a. The corresponding Z sat as a function of volume 
fraction r] = 4irpa 3 /3 is displayed in Fig. 2, with a com- 
parison to its counterpart deduced from the numerical 
solution of PB theory supplemented with the popular 
procedure proposed by Alexander et al [5]. On this fig- 
ure, we have also plotted the results obtained without 
added salt [where the term sinh<£> on the r.h.s. of Eq. 
(1) is replaced by exp<f>, due to the absence of co-ions]. 
Our results are fully compatible with those obtained from 
Alexander's method, with a similar agreement for cylin- 
drical macro-ions (not shown). It is eventually instruc- 
tive to note that for a charged plane confined without 
added salt in a WS slab of width 2h, Alexander's satu- 
ration surface charge may be computed analytically with 
the result [13] cr S at = 2~ 3 / 2 smh(n/V2)/(£ B h) whereas 
we obtain cr sat = \/6Argcosh(5)/(7r^B/i) following our 
prescription. Both expressions agree within 10% and re- 
markably exhibit the same functional dependence on Ib 
and density (through h). 

That our prescription compares favorably with Alexan- 
der's procedure for the planar, cylindrical and spherical 
geometries, calls for the more stringent test of confronting 
our predictions against experimental data. We shall con- 
sider two specific situations corresponding to two differ- 
ent geometries: crystallization of charged colloidal sus- 
pensions and osmotic pressure in B-DNA solutions. 

Crystallization of charged spheres. Investigation 
into the phase diagram of charged polystyrene colloids 
has been conducted experimentally by Monovoukas and 
Gast [8] and compared to the phase diagram of a sys- 
tem where particles interact through a Yukawa poten- 
tial (deduced from extensive molecular dynamics simu- 
lations [16]). However this comparison requires an ad- 
hoc choice for Z c ff. We use here the results we found 
for the effective charge as a function of ionic strength, 
which we insert into the numerical generic phase dia- 
gram of Yukawa systems [16]. We emphasize that there 
is no adjustable parameter in our equations since the ra- 
dius of the polystyrene beads, the only parameter enter- 
ing our description, was independently measured to be 
a = 667 A. We only make the (reasonable) assumption 



that the bare charge Z of the colloids is large enough to 
have Z e ff ~ Z sat . The results for the melting line using 
our prescription for the effective charge are confronted 
to the experimental data in Fig. 3. We also plot the 
result for the melting line for an ad-hoc constant effec- 
tive charge, Z c g = 880, as was proposed in [8] (while in 
our case the latter varies between 500 and 2000 on the 
melting curve). The observed agreement of our results il- 
lustrates both the pertinence of our prescription for Z sat 
and the relevance of the PB saturation picture for macro- 
ions of large bare charge, for monovalent micro-ions in 
water. 
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FIG. 3. Liquid-solid transition of charged polystyrene col- 
loids: volume fraction for melting rj m as a function of salt ionic 
strength 7o. Dots are experimental points for the melting line 
extracted from Ref. [8]. The solid line is the theoretical pre- 
diction for the melting transition using our prescription for 
effective charges. The dashed line corresponds to Z e g — 880 
(see text) . 



Osmotic coefficient of B-DNA. A similar test of 
our method may be performed for the cylindrical geom- 
etry using the experimental results for rigid cylindrical 
polyelectrolytes such as B-DNA [17]. Wc: specifically con- 
sider the measurements of the osmotic coefficient <fr — 
n osm /IT c , defined as the ratio between the osmotic pres- 
sure II osm to the pressure II C of releasable counter-ions 
having bare density c c (H c = ksTcc). Within the WS 
model, B-DNA macro-ions are confined into cylindrical 
cells, whose radius Rws is related to the bare concentra- 
tion of DNA counter-ions as c c = {Idna^Rws)^ 1 ■> w hh 
f-DNA = 1.7 A the distance between charges along DNA. 
The osmotic pressure is related to the densities of micro- 
ions at the cell boundary: n osm = k B T(p + + p_ — 2I ) 
[18], which can be recast in the form II osm = /cbT(kJ — 
K 2 )/(4:Tr£ B ) introducing the screening factor defined 
previously. This latter quantity is computed from our 
prescription, following the same lines as for the spherical 
case [see Eq. (3)]. 
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FIG. 4. Osmotic coefficient of B-DNA vs. density of 
DNA phosphate ions c c , for 7 = 10 mM, 2mM and mM 
(from bottom to top). Dots: experiments of Refs. [17]; solid 
lines: present prescription; dashed line: prediction of Oo- 
sawa-Manning condensation theory. 



lent micro-ions in water at room temperature, micro-ions 
correlations are negligible for all known macro-ions. This 
may no longer be the case in presence of multi-valent 
micro-ions. More generally, PB is a reasonable approxi- 
mation when the macro-ion size a is much larger than £b 
[19], and the saturation plateau of Z e g as a function of 
the bare charge Z is an important physical phenomenon 
that our approach allows to capture. When a becomes 
of the same order as £b the amount of counter-ion "con- 
densation" found in molecular dynamics or Monte Carlo 
simulations is larger than predicted by PB [19,20]; our 
method for Z sat therefore provides an upper bound for 
the effective charge. It is moreover noteworthy that for 
B-DNA where ci/£b — 1.4, our approach still gives a 
valuable first approximation, and that omission of charge 
rcnormalization leads to spurious results (such as nega- 
tive osmotic coefficients corresponding to an unphysical 
phase transition at physiological salt concentrations). 



In Fig. 4, the corresponding results for the osmotic 
coefficient are confronted against the experimental data 
of Refs. [17], showing again a good quantitative agree- 
ment. As in Ref. [18], we report the prediction of classi- 
cal Oosawa-Manning condensation theory, for which the 
osmotic coefficient is constant [cf) — Idna/{^b)] at com- 
plete variance with the experiments. Again we emphasize 
that the only quantity introduced in our description is the 
diameter (a — 10 A) of DNA, known from independent 
measurements. 

In conclusion, we have put forward a simple method 
of asymptotic matching to compute the effective charge 
at saturation for isolated macro-ions. In the situation 
of finite density, this method has been translated into 
a prescription, the validity of which has been assessed. 
This approach (mostly suited to describe the colloidal 
limit na ^> 1) amounts to consider the highly charged 
macro-ions as objects with constant electrostatic po- 
tential <~ AkT/e, independently of shape and physico- 
chemical parameters (size, added 1:1 electrolyte. . . ). As 
a general result we find that the effective charge is an 
increasing function of n, which stems from the reduc- 
tion of the attraction between the counter-ions and the 
macro-ion. Addition of salt consequently brings two an- 
tagonist effects on the effective Coulombic interaction be- 
tween macro-ions: the range of the interaction decreases 
due to screening, while the amplitude increases due to 
the effective charge. The competition between these two 
effects might be a key point in the understanding of these 
systems. 

An important question concerning our approach, is 
that of the validity of PB theory for highly charged 
macro-ions. Within PB, micro-ions correlations are ne- 
glected, but the approach may still allow to describe high 
macro/micro-ions couplings: in particular, for monova- 
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